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Abstract. We investigate the potentiality of using strong lensing clusters to constrain the cosmological parameters 
Qrn and JIa- The existence of a multiple image system with known redshift allows, for a given {Qm, cosmol- 
ogy, absolute calibration of the total mass deduced from lens modelling. Recent Hubble Space Telescope (HST) 
observations of galaxy clusters reveal a large number of multiple images, which are predicted to be at different 
redshifts. If it is possible to measure spectroscopically the redshifts of many multiple images then one can in prin- 
ciple constrain (f2m, ^a) through ratios of angular diameter distances, independently of any external assumptions. 
For a regular/relaxed cluster observed by HST with 3 multiple image systems, each with different spectroscopic 
redshifts, we show by analytic calculation that the following uncertainties can be expected: Qm = 0.30 ± 0.11, 
= 0.70 ±0.23 or Qm = 1.00 ±0.17, ^a = 0.00 ±0.48 for the two most popular world models. Numerical tests 
on simulated data confirm these good constraints, even in the case of more realistic cluster potentials, such as 
bimodal clusters, or when including perturbations by galaxies. To investigate the sensitivity of the method to dif- 
ferent mass profiles, we also use an analytic "pseudo-elliptical" Navarro, Frenk & White profile in the simulations. 
These constraints can be improved if more than 3 multiple images with spectroscopic redshifts are observed, or by 
combining the results from different clusters. Some prospects on the determination of the cosmological parameters 
with gravitational lensing are given. 
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1. Introduction 



A new "sLandaid cusiuulugiual uiudd" lias aiisun lii Lliu 
last few years, favoring a flat Universe with Qm ^ 0.3 
and il\ ^ 0.7. This is mainly based on the results of two 
experiments which give roughly orthogonal constraints in 
the {^lrn,^\) plane (see Fig. ^ for a recent update). The 
first one is obtained by considering type la supernovae 
(SNIa) as standard candles. The detection of a sample of 
high redshift SNIa (up to z ~ 1) by two groups favours 



a non- vanishing cosmological constant (Pcrlmutter et al 



199S; Ricss et al. 1998), large enough to produce an accel- 
eratin g cjcpanoion. However, ovidonoo for a non aero coo 



favor a flat Universe ( Balbi et al. 2000 ; Melchiorri et al 
2000 ). However, there still remains a degeneracy in the 
combination of flm and Qx because CMB experiments are 
primarily sensitive to the total curvature of the Universe. 
Even with the accuracy of the future MAP and Planck 
missions, the constraint issued from the CMB alone will 
be degenerate. 

The combination of these two sets of constraints has 
led to the currently favored model of low matter den- 
sity and a non-zero cosmological constant, preserving a 
flat geometry 



mologi bal constant ia atill oontrovoroial, oinoo aupornovao 
might evolve with redshift and/or may be dimmed by in- 
tergalactic dust ( Aguirre 1999D . The fundamental assump- 
tion of a homogeneous Universe and its implication for a 



g- 


White 199?; 


Efstathiou et al. 1999; 


Sahni & Starobinsky 200C; Jaffe et al. 



non-zero cosmological constant are also discussed ( Celerier 



200C 



Kolatt & Lahav 2001). The second constraint is 



derived from the location of features in the cosmic mi- 
crowave background (CMB) anisotropy spectrum, partic- 
ularly the flrst Doppler peak. The most recent results ob- 
tained with the Boomerang and MAXIMA experiments 



2001). Although these recent results are quite spectacular, 
there still remain many sources of uncertainties with both 
methods. Thus any other independent test to constrain 
the large scale geometry of the Universe is important to 
investigate. Gravitational lensing, an effect involving large 
distance scales, has been considered as a very promising 
tool for such determinations. Indeed, the statistics of grav- 
itational lenses depend on the cosmological parameters via 
angular size distances and the comoving spherical volume 
(e.g. [Turner et al. 198^ ; [Turner 1990| ; [Kochanek 1996 ; 
Falco et al. 199^ ). This technique has provided an upper 
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1999), lensed radio sources ( Cooray 1999), lens ed galaxies 
in the Hubble Deep Field ( [Cooray et al. 1999| ). Although 
most authors favor a lambda-dominated flat Universe, 
there remain some uncertainties in the mass distribution 
of the galaxy lenses and on the luminosity function of the 
sources. Evolutionary effects may also play a role in these 
statistics. 

Another application of gravitational lensing to con- 
strain the cosmological parameters is to use the statis- 
tics of the "cosmic" shear variance. Van Waerbeke et al 



(1999|)]showed that it is related to the power spectrum of 



the large scale mass fluctuations, and then to fin 



,. The first 
in the 



rang 


e 


0.2-0.51 ( 


Maoh et al. 2001 


7 


Van Waerbeke et al. 


2001 


; 1 


Jacon et al. 2000; 


Wittman et al. 2000 


). Imaging 



surveys with the next generation of panoramic CCD cam- 
eras will reinforce this very promising technique. In the 



case of weak lensing by clusters of galaxies, Lombard! & 



Bertin (1999 ) and Gautret et al. (2000 ) suggested methods 
to constrain the geometry of the Universe. These methods 
need however to recover the mass distribution and/or to 
know acurately the redshift of a huge number of distant 
galaxies, making this method not practical in the near 
future. 

In this paper, we focus on a measurement technique of 
[firm fl\) using gravitational lensing as a purely geomet- 
rical test of the curvature of the Universe, since the lens 
equation depends on the ratio of angular size distances 
which is sensitive to the cosmological parameters. In the 
most favorable case, a massive cluster of galaxies can lens 
several background galaxies, splitting the images into sev- 
eral families of multiple images. The existence of one fam- 
ily of multiple images, at known redshift, allows to cali- 
brate the total cluster mass in an absolute way. In the case 
of several sets of multiple images with known redshifts, it 
is possible in principle to constrain the geo metry of the 
Universe. This method was pointed out by 



Blandford & 



Naray ou (1992 ), and earlier suggested by [Paczynski fc 
Gorski (1981 ), but the uncertainties in any lens studies 
were considered too large compared to the small varia- 
tions i nduced by the cosm ological parameters. More re- 
cently. Link fc Pierce (199 j ) (hereafter LP98) re-analysed 
the question in the light of the typical accuracy reached 
with HST images of clusters of galaxies. Following their 
method, which inspired our work, we try to quantify in 
this paper what can be reasonably obtained on (fljmflx) 
from accurate lens modeling of realistic cluster-lenses. 

The paper is organized as follows. In Sect. 2 we sum- 
marize the main lensing equations and we introduce the 
relevant angular size distance ratio which contains the de- 
pendence on the cosmological parameters. The variation 
of this ratio is then compared to that of other variables 
(lens potential parameters and redshifts) to derive the ex- 
pected uncertainties on ri,„ and fl\. In Sect. 3 we present 
the method in detail and the results from simulations of 
various types of families of images and of different types of 
lens potentials. Some conclusions and prospects for the ap- 
phcation to real data are discussed in Sect. 4. Throughout 




0.8 1.0 



Fig. 1. Constraints on {Vim,fl-\) derived from the most 
recent results of the BOOMERANG and MAXIMA-I ex- 
periments on the CMB fluctuations and the last results 
from the SNIa analysis (from Jaffe, 2001). Overplotted is 
the ratio of the geometrical factor of the lens equation 
for two source redshifts, 6(251,252) = E{zsi) / E{zs2) as 
discussed in the text. In this example, the redshifts are 
chosen as representative of a typical lens configuration: 
zl = 0.3, 2:51 — 0.7 and 252 — 2. 



this paper we assume Hq = 65 km s~^ Mpc~^ (note how- 
ever that the proposed method and results are indepen- 
dent of iJo). 

2. Influence of Q,m and Q,x on image formation 

2.1. Cosmology dependent lensing relations 

We first define the basic mathematical framework, follow- 
ing the formalism presented in Schneider ct al. (1992| ). We 
consider a lens at a redshift zl with a two-dimensional 
projected mass distribution S(0) and a projected gravita- 
tional potential (f>{9), where is a two-dimensional vector 
representing the angular position. A source galaxy with 
redshift 23 is located at position 6s in the absence of a 
lens, and its image is at position Oj. In the lens equation 



es = Oi w^{ei) 



(1) 



DoL, Dls and Dos are respectively the angular diameter 
distances from the Observer to the Lens, from the Lens to 



the Source and from the Observer to the Source (Peebleg 
1993). if is the reduced gravitational potential which sat- 
isfies V^(p = 2 E/Scrit with the critical density 



S 



Dos 



crit 



4nG DlsDol 



(2) 
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In these equations the dependence on the cosmologi- 
cal parameters appears only through the angular diameter 

{E as "e/- 



distances ratios F — ^'^^ and E 



Dos Dos 
ficiency") for a fixed cluster redshift. They correspond to 

a scaling of the lens equation, reflecting the geometrical 

properties of the Universe. 

In the general case, we can scale the potential gradient 

as: 



V0(Qj) -up^ DoL f{Oi,0c,a,. 



Note that the discrepancy between the different cos- 
mological parameters is not very large, less than 3% be- 
tween the Einstein-de Sitter model (EdS) and a flat, low 
matter density one. Moreover, a characteristic degeneracy 
appears in the {Qm, ^x) plane, which is roughly orthogo- 
nal to the one given by the detection of high redshift su- 
pernovae, and quite different from the CMB constraints. 
A similar degeneracy was also found in the analog weak 
lensing analyses, (Lombardi & Bertin 1999) or by Gautret 
et al. (200q). 



where (Tq is the central velocity dispersion and / a di- 
mensionless function that describes the mass distribution 
of the cluster. It can be represented by flducial parameters 
such as a core radius 6c or a mass proflle gradient a. The 
lens equation reads 



0s — Oj 



4 Dis 

C2 



OS 



f{ei,ec,a,...) 



(4) 



Oi-^ fiei,9c,a, . . . ) X Ein„,, fix, ZL, zg) 



We will focus on the E'-term which entirely contains 
the dependence on Qm and fl\ and which is independent 
of Ho. 

2.2. The E-term 

For a given lens plane zl, the ratio E increases rapidly 
as the source redshift increases, and then flattens at zg ~ 
1.5. There are also small but significant changes with the 
cosmological parameters (Fig. ||). The dependence of E 
with respect to flm is weak, whereas the variation with 
fl\ is larger. 

We now consider fixed redshifts for the lens and 
sources. Assuming a fixed world model, a single family of 
multiple images can in principle constrain the total cluster 
mass as well as the shape of the potential, removing the 
unknown position of the source 6s using Eq.(^). In prac- 
tice, good constraints on the shape of the potential / are 
obtained with triple, quadruple or quintuple image sys- 
tems. However the absolute normalization ctq of the mass 
is degenerate with the E-term, that is with respect to Q„i 
and fix. 

2.3. Ratio of E-terms for 2 sets of source redshifts 

To break this degeneracy a second family of multiple im- 
ages is needed. To get rid of the strong dependence on ctq, 
it is useful to consider the ratio of the positions: 



E{zs{)\\f{ej.,...)-f{ej.,...)\\ 

\Oil - E{zs2)\\f{ei.,. ..)- . . . )|| 



(5) 



(here and hereafter, the superscript refers to a family and 
the subscript to a particular image within a family) . This 
ratio is plotted in Fig. |lj highlighting the infiuence of flm 
and fix through the ratio 6(2:51, Z52) — E{zsi) / E{zs2)- 



Another approach to quantify the dependence of a 
given lens configuration on flm and fix is to fix the lens 
redshift and to search for two source redshifts zsi and 
Z52 which give the largest variation of the i?-term when 
scanning the {flm,^\) plane. For illustration we arbi- 
trarily choose two sets of cosmological parameters, for 
which the relative variation of E is large: CPl {flm = 
0.3, = 0) and CP2 {fl^ = l^fix = 0, i.e. the EdS 
model). Varying zsi and zs2i the function e{zsi,zs2) = 
ecP2{zsi, ZS2) / ecpi{zsi, ZS2) - 1 represents the percent- 
age of discrepancy between CP\ and CP2 for zsi > 
ZS2 (> ^l) (Fig. I). 

For a given high-redshift zs2 source the best lowest 
source redshift is zsi — zl, and for a given zsi the best 
zs2 is the highest redshift, the difference between cosmo- 
logical models increasing with zs2- In all cases, this rela- 
tive difference is of the order of a few %, meaning that the 
lens mass distribution must be known to the same degree 
of accuracy to get further constraints on the cosmological 
parameters. Hence, for 2 systems of images, the best con- 
figuration is one background source close to the lens, in 
the rising part of E{zs) and another one at high redshift, 
to take into account the asymptotic value of the ratio. 
Note however that for a source redshift close to the lens, 
the i?-term becomes very small. Also, the location of the 
images is very close to the lens center which makes the 
detection of multiple images quite improbable, as small 
caustic sizes imply small cross sections. 

2.4. Relative influence of the lens parameters 
2.4.1. Physical assumptions 

In order to quantify the expected uncertainties on flm and 
fix, it is possible to analytically estimate the influence 
of the different lens models parameters. We use a model 
of the potential derived from the mass density described 



by Hjorth & Kneib (2001), hereafter HK. It is based on 
a physical scenario of violent relaxation in galaxies, also 
valid in clusters of galaxies. The mass density is charac- 
terized by a core radius a and a cut-off radius s: 



p{r) 



Po 



1 



(6) 



Then the projected mass density S(6') is represented by: 



no) - So 



(7) 
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Fig. 2. Variation of the lens efficiency -^(^^s) as a function of the source redshift for different sets of cosmological 
parameters:. (Left): JIa = and Vim varies from 0.1 to 1.0. (Right): Vim = 0.3 and Vl\ varies from to 0.7. In both 
cases, the lens redshift is zl = 0.3. For clarity, above each main plot, we also plotted the same curves, normalised with 
the E{zs) function for {Q.m,Vlx) = (0.3,0.). 
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Fig. 3. Relative difference between the ratio e(zsi, Z52) = 
E{zsi) / E{zs2) for two extreme cosmological models: 
{Vim = 0.3, = 0) and {Vtm = l,rjA = 0). This func- 
tion e{zsi, ZS2) is plotted in the {zsi, ZS2) plane, assuming 
zsi < zs2 and zl = 0.3. 

where 9 is the angular coordinate, 9a = a/DoL and 9s — 
s/ DoL- So is a normalization factor, related to the cluster 
parameters: 



irpoD 



OL 



(8) 



Finally,the mass inside the projected angular radius 9 is: 

fs — fa 

^9^+91-^/9^+91 + 9s -9a\. (9) 

The velocity dispersion (t{9) is related to the mass den- 
sity and the gravitational potential via the Jeans equation. 
Assuming an isotropic velocity dispersion and retaining 



terms up to first order in 9a /9s, we get the relation be- 
tween the central velocity dispersion ctq = c(0) and pq: 



pos. 



(10) 



Finally we compute the expression of the deviation angle 
between the positions of the source and of the image due 
to the lens: Dg^ — \\9i — 9s\\, neglecting second order 
terms in 9a /9s (we suppose here s ^ a): 



Ds, = 



16 ^2 



1 r 



V9^+9l - V9^+9l + 9s - 9a 



(11) 



2.4.2. The single multiple-image configuration 



The central velocity dispersion (or equivalently the mass 
normalization of the cluster core) is obviously the pre- 
dominant factor in any lens configuration. With a sin- 
gle family of images we can only constrain the combi- 
nation cr§i? and cannot disentangle the influence of the 
cosmological parameters and the absolute normalization 
of the mass (Fig. If we were able to measure the to- 
tal mass within the Einstein radius independently from 
lensing techniques and with an accuracy better than a 
few %, we could in principle put some constraints on f2A. 
Observationally, there are 2 situations where it is likely 
that we could disentangle the effect of cosmology and ab- 
solute mass: 

1) in the case of a cluster- lens with extremely good X-ray 
data, particularly in estimating the temperature distri- 
bution of the gas (under the assumption of hydrostatic 
equilibrium) , 

2) in the case of a multiple system around a single galaxy, 
for which one is able to measure accurately the stellar 



velocity dispersion of the lensing galaxy (Tonry & Franx 
1999|). 
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Fig. 4. Variation of the central velocity dispersion cto in 
the {Q.m,^\) plane, assuming that the product aQE{zs) 
is constant, uo has been fixed to 1000 km s^^ for Vim — 1, 
f^A = while zl = 0.3 and zs = 1. 

However in both cases this represents some observational 
challenge and requires the most powerful instruments to 
achieve this goal. 

Although the error budget in the image positions is 
dominated by the error on the total cluster mass (or equi- 
valently the velocity dispersion), we can determine the 
relative influence of the other parameters to infer the im- 
portance of rim and Vl\ in the image formation. The rel- 
ative error on the deviation angle De^ depends on ctqj (^a, 
0s and 9 1 for the gravitational potential, and zl, -^Sj 
and fix for the i?-terni (Eq.(|ll|)). Therefore we can write: 



dE 



d0„ 



d0T 



a.01 ■ 



(12) 



with 

d£: 



dVlr, 



dVl\ dzL , dzs 
OiQ.^-^ ha^L l-"zs • 



(13) 



ag^ , ae^ and ag^ can be computed analytically while 
= 2 is the largest factor. Since the angular diame- 
ter distances do not have an analytic expression if fix is 
non-zero, the coefficients ao„, ct^ixi ctzi, and a^g must be 
computed numerically. In practise, they are computed for 
a given set of parameters (zi^,Vlm,Vl\) as their variation 
with Vim and Vl\ is of higher order. For a reasonable set 
of lens parameters, the a-coefficients are of the same or- 
der of magnitude, except that a^g and can dominate 
the error budget if the source redshift is close to the lens 
(an unlikely case). On the contrary, an^^ and aux are of 
second order, and a^^ is somewhat larger than ao^. This 
reflects again the fact that i?-term is more sensitive to Vl\ 
than to Vim- 

To quantify the relative influence of all the parameters 
in the case of a single family of images, we computed ex- 



plicitely dDgjDgj in two cases, for a cluster- lens and for 
a galaxy- lens. 

1) For a cluster of galaxies, we take the following param- 
eters: ZL = 0.3, zg = 4, Vim = 0.3, Vlx = 0.7, Os/ea = 10 
and Oi/Oa = 4. We thus find from Eq.(|l): 



dDg 



2^-0.31^ 



0.21 



dOn 



Dg, (To 

- 0.17— + 0.062— + 0.012 ^ 

zl 



-0.51 

dVlr. 



dOs 



0.14 



dVly 



Vlx 
(14) 



Let us assume a perfectly known mass profile (i.e. 
d9a = d9s = 0.). Neglecting the influence of Vim, we ask 
what precision would be required on ctq to derive an error 
of 50% on Vlx- The accuracy of the position of the center of 
the images is calculated using the flrst moment of the flux 
fi9) on a given image: 0i = J 9f{0)d9^ / J f{9)d9'^ which 
yields an error d9i of a fraction of the spatial resolution. 
HST observations are then required to reach d9i = 0.1" 
(LP98) or better. To reduce the uncertainty on the red- 
shift measurements, we assume spectroscopic determina- 
tions, so that dz ~ 0.001. Finally we have to compute the 
relative errors on Dgj = — so the position of the 
source is in principle required. But as we are in the strong 
lensing regime, we assume that 9s <C 9i, so that both 

dDg. d9i 

quantities Dg^ ~ W/ and — — ~ — — are directly related 

Dgj 9i 

to observable ones. Taking these values into account, we 
need to know ctq with 3.6% accuracy to get the expected 
constraint on Vlx- Such an accuracy is out of reach with 
observations of clusters of galaxies. 

2) For a single galaxy, we consider typically: zl = 0.3, 
zs = 4, Vim = 0.3, Vlx = 0.7, 9sl9a = 200 and 9il9a = 
200 (ratios take n from the modelisatio n of the lens HST 
14176+5226 by [Hjorth fc Kneib (200l| )), leading to: 



dPg, 
Dg, 



2^-0.72^- 0.0066^ + 0.73^ 

9l 9a 9s 



0.17^+0.062^ + 0.012^ 

Zl zs i2m 



■0.14 



dVlx 

nx 

(15) 



Taking the same values for the observational errors 
and considering a perfectly known mass profile, we require 
an accuracy of 6.4% on ctq to derive a 50% error on Vlx- 
For a typical galax y, this represents about 15 km s~^. 
Warren et al. (199^ ) measured the velocity dispersion in 



the deflector of the Einstein ring 0047-2808 with an error 
of 30 km s~^. A better accuracy could be obtained by 
looking at particular strong absorption features with 10m 
class telescope observations. This could be sufficient to 
confirm an accelerating Universe. 

2.4.3. Configuration with 2 multiple-image systems 

With a second system of multiple images another region 
of the E(zs) curve is probed while the cluster parameters 
are the same. In that case, the relevant quantity becomes 
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the ratio of the deviation angles for 2 images 6ji and 6j2 
belonging to 2 different families at redshifts 2:51 and 232- 

We define Rg , e — — '—■ This function has the advan- 
De^, 

tage that it does not depend on uq anymore. Following 
our previous definitions, we can write 



di?« 



<iE{zs2) dE{zsi) 



Re, 



E{zs2) 
aej {O12) 



E{zsi) 

d6l/2 



- ae, {Oil ) 



d0. 



in 



+ (aeAOi^)~aeAOn)) ^ 



7/2) 



^71 



)) 



(16) 



Numerically, we chose a typical configuration to com- 

pute ejOa - 10, Oi2/0a = 4, 012/eii - 2, 

^6^1,61^2 

zl = 0.3, zsi — 0.6, ZS2 — 5, assuming fl„i = 0.3 and 
f^A = 0.7. This gives the following error budget: 



= 0.92- 



dzT 



0.99- 



dz, 



SI 



0.062 



0.018 



ZL 

d0ii 



- 0.31 



ZSI 

d0j2 



0.12- 



dzs2 

ZS2 

d0„ 



7/1 



7/2 



0.21. 



d0.. 



0.034 



dQr 



0.037 



djh 



(17) 



The contribution of the physical lens parameters in 
this error budget is strongly attenuated comprared to the 
single family case. There is no more dependence on cto and 
the dependence on the mass profile ( 0a, Os) is reduced by 
about a factor of 2 compared to a single family of images. 
This corresponds to the variation of the potential between 
011 and 012, the absolute normalization being removed. 
Anyhow, this can still represent the main source of error 
because we cannot expect to constrain 0a to better than 
1.5% and 0s to better than 2% typically (see Sect. |3.3.3|) . 

For the source redshifts, we have selected one of the 
sources at zsi — 0.6, which means that its a-coefficient is 
quite large. The accurate value of the redshifts is thus fun- 
damental, and a spectroscopic determination is essential 
{dz ~ 0.001). A photometric redshift estimate would not 
be satisfactory, because we cannot expect an accuracy bet- 
ter than 10% in most cases (dz : 

(2000|)]. We keep d0i = 0.1". The strong lensing regime 

. , , ^ ||6»j2-6»s2 
approximation leads to Re^i ,9^2 

d^0ji,e^2 d0j2 d0ji 



\e, 



6»c 



Uj2 

- — and 



Re, 



012 



We can then separate the contributions of the param- 
eters that do not depend on Q,n or fl\ from those which 
depend on them and re- write Eg. (p^) : 



Ac 



A, 



dflx 



(18) 



'Errl^ {0J^ ,0j2,0a,Os)+ Err22 (O^, f]^, ^L, ^51,^52) 



An„ and Aq^ depend on fim, ^^a, -^l, zsi,zs2 while Errl^ 
and Err2^ are the quadratic sums of the errors, with a 
separation between the geometrical parameters and those 
depending on the cosmology. For each set of cosmological 
parameters we then compute all these coefficients numeri- 
cally. In addition, we also need a calculation of the "degen- 
eracy" dilm/dilx to obtain either dilm or d^l\. This is the 
slope of the degeneracy curves of the E-terms ratio plotted 
in Fig. ^ Indeed considering 2 points (11m, ^^a) and {ilm + 
dflrm ^\-\-dfl\) on such a curve (for a given set of zl, zsi, 
Z52), we have e{il.m,^\) = e(il„i + dfl„i,il\ + dfl\) so 
that we get dn,n/dQ.\ = -dn^e{zsi, zs2)/dn,^e{zsi, zs2)- 
The final expected errors on dQm and dQ\ are plotted in 
Fig. H for a continuous set of world models. The method 
is in general far more sensitive to the matter density than 
to the cosmological constant, for which the error bars are 
larger. This apparent contradiction with the general state- 
ment that lensing is more sensitive to the cosmological 
constant than to the matter density is due to the fact that 
we analysed the ratio of two E-terms and this ratio varies 
more rapidly with when scanning the {^m, ^\) plane 
(Fig|^). For illustration, we quantitatively obtain the fol- 
lowing errors for the corresponding cosmological models: 



ACDM 
SCDM 



Sn„i =0.11 
SQr,, = 0.17 



Silx = 0.23 
Snx = 0.48 



This analysis shows that the expected results are quite 
encouraging, and the constraints we could get are similar 
to the ones currently obtained by other methods. Note 
however that these typical values require both HST imag- 
ing of cluster lenses and deep spectroscopic observations 
for the redshift determination of multiple arcs. They may 
depend on the choice of the lens parameters and on the 
potential model chosen to describe the lens, a problem 
that we will now investigate. 

3. Constraints on the cosmological parameters 
from strong lensing 

3.1. Existence of multiple systems of lensed images 

More and more cluster-lenses are known to show several 
systems of multiple images (with spectroscopic or photo- 
metric redshifts). Lens modelling is then performed with 
a good accuracy and allows the prediction of extra fam- 
ilies of images and their expected redshifts. When these 
images are later identified and if their redshift can be mea- 
sured spectroscopically, an iterative process brings the lens 
model to a high level of accuracy, where most of the pa- 
rameters which characterise the mass profile are strongly 
constrained. This full process has been applied success- 



fully in a few clusters such as A2218 (Kneib et al. 1996) 



A370 (Kneib et al. 1993a; 3mail et al. 1996; Bezecourt 



3t al. 1999a) or AC114 (Natarajanet al. 199S; Campusanc 
3t al. 200 1| ). 

In order to apply more systematically the method pro- 
posed here, we may ask whether the few cited cluster- 
lenses are representative of some generic cluster or if they 




correspond to very peculiar configurations. To answer this 
question, we simulated a typical cluster at redshift z — 0.2 
with the following characteristics. A main clump is de- 
scribed with the potential of Eq.(^), the so-called HK mass 
density, with a = 50 kpc and s — 500 kpc. These values 



are typical of cluster-lenses at this redshift ( pmith et al. 



200 1|). |The central velocity dispersion is varied from 800 
to 1400 km s^^ to allow a variation of the Einstein radius. 
In addition, 12 individual galaxies are added in the mass 
distribution, following the prescription used by Natarajan 
& Kne b (1996 ) and for a total contribution of 30% of the 
total mass. Their individual masses are scaled with respect 
to their luminosity Li by the Faber-Jackson relation: 



(19) 



with =150 km s ^ (following Faber et al. (1997 )), and 
with a cut-off radius: 



(20) 



providing a constant ratio MjL ( Natarajan fc Kneib 
19970 

To simulate the background sources, we used the 
Hubble Deep Field (HDF) image acquired by the ES T 
( IWilhams et al. 1996]) . From |Fernandez-Soto et al. (1999| ) , 



946 galaxies were extracted from the deepest zone of the 
F814W image, up to a magnitude limit ^^(8140) = 28.0 
and over an angular area of 3.92 arcmin^. These authors 
also provide a catalog of photometric redshifts for all these 
objects. In addition, for about 10% of them, a spectro- 
scopic redshift is available. We used this redshift distri- 
bution (spectroscopic redshift preferably used when avail- 
able) as a sample of galaxy-sources to be lensed by the 
simulated clusters. In order to increase the statistical sig- 
nificance of this simulation, we generated a source cat- 
alogue with 10 times the number of galaxies extracted 
from the HDF image. We then distributed these sources 



at random angular positions over the central inner 40 x 40 
arcsec^. We checked that this region includes the external 
radial caustic line, so that no multiple images are lost. The 
increase in the galaxy density is then corrected for in the 
final results. 

Table |l] presents for each value of the central velocity 
dispersion the number of systems found with their image 
multiplicity. We also determined the number of systems 
in which each image could be observed (with a magni- 
tude ^i?(8140) < 24.5, corresponding to typical HST in- 
tegration time of 10 ksec). Objects with ^^(8140) > 28.0 
could be observed due to the lens effect if the magnifica- 
tion exceeds a factor of 25. This very rare configuration 
is neglected in our simulations for simplicity. For a clus- 
ter massive enough (co > 1200 km s^^, corresponding to 
^'^tot > 2. 10^"^ Mq for our potential model), numerous sys- 
tems of multiple images (mainly triple images) are formed 
and a significant fraction could be observable. Although 
these simulations are quite simple and cannot be used for 
realistic statistics of image formation, it gives us confi- 
dence that the use of multiple image families for the de- 
termination of the cosmological parameters is achievable 
and should be applied on a large number of rich clusters. 



3.2. Method and algorithm for numerical simulations 

In most cases, clusters of galaxies present a global elliptic- 
ity in their light distribution or in their gas distribution 
traced by X-ray isophotes. It is generally believed that this 
is related to an ellipticity in the mass distribution. This 
has indeed been recognized several times by the modeling 
of cluster lense s such as MS2137 -23 dMelUer et al. 1993D 
or Abell 2218 (Kneib et al. 1996). So we include such an 
ellipticity in our modeling of cluster potentials. The basic 
distribution of matter we consider is again the HK one, 
with, in addition, a substitution of the radial distance r 
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Table 1. Number of systems of images obtained for sim- 
ulated cluster potentials with different values of the cen- 
tral velocity dispersion ctq (the corresponding Einstein ra- 
dius Re is given for z = 1). The redshift distribution 
of the sources is assumed from the HDF data. We take 
random positions for the sources over the central inner 
40 X 40 arcsec^. nj is the number of systems of j images. 
n* is the number of systems of j images with ^5(8140) < 
24.5, corresponding to "observable" ones. Then, each sys- 
tem is counted both in rii and in n* with j < i (if only 
j images among i are detectable). Systems counted in ng 
show no "observable" image. So ni-|-ntot = + '^■1 + ?^tot ' 
which is the number of galaxies in the selected field. 





(km s"^) 


800 


1000 


1200 


1400 


Re 


(arcsec) 


5 


14 


28 


40 



"■0 

* 

no 



(1) 



78 




73 




69 




65 




(2) 

I 

»(3) 



107 
29 



107 

34 



99 
34 



66 
26 



-0 — 0X168- 



O.M 



-OlO- 



-il3- 



_QJJ QAL 



n3 



0.12 
0.034 



0.60 
0.068 



8.0 
1.7 



41 
3.9 



714 0.034 
nl 0.011 



0.011 
0.011 



0.034 
0.034 



0.057 
0.011 



nt 



0.022 
0.011 



0.011 




0.11 





0.011 
0.011 



ne 0.011 

ng 



0.011 
0.011 



7 


717 0.011 


0.011 


0.011 






0.011 








8 


na 


0.011 










* 

rig 











total (j > 2) 


Tttot 


0.18 0.70 


8.3 


41 






0.17 0.50 


4.0 


17 


including 0.6 


galaxies at 


z < 0.2. 






'■^Hncluding 5.2 


galaxies at 


z < 0.2. 






including 4.6 


galaxies at 


z < 0.2. 







by R defined as: 



R 



XcosO + Y sin t 
1+1 



-X sin 6* + y cos ( 



(21) 



where X = [x — Xg) and Y — (y — j/o)- The potential (j) 
is then characterized by 7 parameters, namely: xq, yo, e, 8 
for the geometry of the lens and aQ,9a, ds for the shape of 
the mass profile. 

Another popular density profile to be tested is the 
so-called Navarro, Frenk & White (NFW) density pro- 



file found in many simulations of dark matter and cluster 
formation (Navarro ct al. 1997): 



p{r) 



Pc 



(22) 



where pc is a characteristic density and a scale radius. 
No analytic developments have been proposed so far for 
the corresponding ellipsoidal profile. In a companion pa- 
per (Golse & Kncib 2002) we propose a new "pseudo- 
elliptical" NFW profile and compute its lensing properties. 
The corresponding potential is characterized by 6 param- 
eters: xo,yo, e, 6 for the geometry of the lens and Vc, 9s for 
the shape of the mass profile. The characteristic velocity 
Vc is defined by 



(23) 



as explained in Golse & Kneib (2002). 

To create a simulated lens configuration we need to 
fix some arbitrary values of the cosmological parameters 
(ri^,ri°) as well as the cluster lens redshift zl- The nu- 
merical code LENSTOOL developed by one of us (Kncit 



1993D can then trace back the source of a given image or 



determine the images of an elliptical source galaxy at a 
redshift zg. The initial data are several sets of multiple 
images at different redshifts. In all cases we do not take 
into account the central de-magnified images, which are 
generally not detected. With these observables, we can re- 
cover some parameters of the potential while we scan a 
grid in the (17™, fi^) plane. The likelihood of the result 
is obtained via a ^^-minimization (with a parabolic or a 
Monte Carlo method) , where is computed in the source 
plane as: 



X 



EE 

4=1 J = l 



(24) 



The superscript i refers to a given family of multiple im- 
ages and the subscript j to the images inside a family of 
Ui images. There is a total of '^i = images, and 

Sr=i 2('T'j ^ 1) = constraints on the models assuming 
that only the position of the images are fitted. Qgi is the 
source position associated with the image 0ji in the lens 

3 

inversion. ©sG* is the barycenter of all the Q ai belonging 
to the same family i. Aj is the magnification matrix for a 
particular image and aji is the error on the position of the 
center of image Oji. Quantitatively we will take aj = 0.1" 

3 

for all images, assuming that their positions are measured 
on HST images. 

computed from Eq. (^J) in the source plane is math- 
ematically equivalent to computed in the image plane, 
written as: 



EE' 

i=i j=i 



(25) 
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Table 2. Details on the 3 sets of multiple images used in 
the simulations in Sect. 3.3.1 and 3.3.2. Ui represents the 



number of images used for each family. It does not include 
the central de-magnified image created for tangential im- 
ages. Nc is the number of constraints in the lens modeling 
for each family. Nc = 2 x — 2 for x and y position. The 
unknown position of the source {xs,ys) is then removed, 
reducing Nc by 2 units. 



Family 


Type 


rii 


2S 


Nc 


i = 1 


Tangential 


4 


0.6 


6 


i = 2 


Radial 


3 


1. 


4 


i = 3 


Tangential 


4 


4. 


6 




where 6j(ji is the image of ©sg* close to 6ji . Indeed Ogi — 
9sc^ = SSj and 6ji — djc^ = ^Ij are assumed to be 
small quantities compared to the variation scale of the 
elements of the magnification matrix A'j. Therefore the 
local transformation from the image plane to the source 
plane is written as <5 Jj = Aj SSj . The main motivation 
for working in the source plane is numerical simplicity 
because the mapping from the source to the image plane 
is not a one-to-one mapping and we may not recover all 
the images when solving the lens equation. 

If Alp is the number of fitted parameters for the poten- 
tial, there is a total of M = Mp + 2 adjustable parameters 
(including and and Nc independant data points. 
We compute a x^-distribution for = Nc — M degrees of 
freedom. In practice, in our simulation we try to recover 
only the most important parameters, like ao (or (Tc), 9a or 
9s, to limit the number of degrees of freedom. This would 
be the case in a real application. 



3.3. Numerical simulations in different configurations 

To recover the most important parameters of the poten- 
tial, we generated 3 families of multiple images (2 tangen- 
tial ones and a radial one for a total number of constraints 
V — 16, see Fig. || and Ta ble |^) with the pseudo -elliptical 
NFW profile developed in |Golse fc Kneib (2002| ). We also 
chose regularly distributed source redshifts (Table ^). The 
4 geometrical parameters of the cluster lens were left fixed 
during the minimization xq = yo ^ 0, 6 = 0° and e = 0.1), 
while the 2 parameters of the potential {vc and 9s) were 
allowed to vary. The initial values for these parameters, 
used to create the set of images, correspond to reasonable 
values found in cluster lenses: ds = 31.3" (i.e. 150 kpc) 
and Vc = 2000 km s~^. This last value corresponds to a 
"classical" central velocity dispersion ctq = 1230 km s~^ 



Fig. 6. Multiple images generated by a pseudo-elliptical 
NFW cluster at zl — 0.3 with the lens parameters: Vc = 
2000 km/s, 9s = 31.3"(rs = 150 kpc) and e = 0.1. Close to 
their respective critical lines, 3 families of multiple images 
are identified: a tangential one (# 1, zgi = 0.6), a radial 
one (# 2, zs2 = 1) and another tangential one (# 3, Z53 = 
4). Units are given in arcseconds. 

3.3.1. Simple cluster potential 

In this case, the number of degrees of freedom is = 
16 — 4 — 12 as 2 cluster parameters are fitted. The con- 
fidence levels of the minimization are plotted in Fig. ^ 
The trajectory of the minimum includes the initial point 
(0.3,0.7) in the (Q^nM^) plane with Xmin = 0. The de- 
generacy in the cosmological parameters is found as ex- 
pected in Fig. |l|. Tighter constraints can be deduced on 
flrn than on fix. We also recover the cluster parameters 



quite satisfactorily with: Vc 



2000;^!^ km 



for a HK model (see Sect. ^.3.21 ). 
the ACDM values (0.3,0.7). 



(f7^,f70) were fixed to 



(Fig. 1) 

and 9a — Note that these errors represent only 

the variations of the fitted parameters when we scan the 
(flrm^^x) plane during the optimisation process. 

This preliminary step corresponds to the "ideal" case 
where we recover the same type of potential we used to 
generate the images. Moreover, the morphology of the 
cluster is regular without substructure, and we included 
one radial system among the families of multiple images. 
These images are known to probe the cluster core effi- 
ciently. Finally, the redshift distribution of the sources is 
wide and the selected redshifts are well separated, for an 
optimal sampling of the E-term. One could ask whether 
any such lens configuration has already been detected 
among the known cluster lenses. It seems that the case of 
MS2137.3-235 3 (;zl = 0-31) is qu ite close to this type of 
configuration (Mellier et al. 1993) with at least 3 families 
of multiple images, including a radial one. Uunfortunately, 
no spectroscopic redshift has been determined for any of 
the images so far. 
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1.0 



Fig. 7. confidence levels in the {Vim,^\) plane ob- 
tained from the optimisation of the lens configuration 
shown in Fig. |^. The 2 main cluster parameters Vc and 
6s were recovered with x^jn = 0. The cross (+) represents 
the original values (I^S^ii^a) = (0.3,0.7). Dark to light 
colors delimit the confidence levels (from 1-a to 4-ct). 
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Fig. 9. confidence levels in the (rim,r2A) plane ob- 
tained from the optimisation of the lens configuration 
shown in Fig. ^. In this plot, the potential was fitted with 
a model different from the initial one (an elliptical HK 
profile instead of a pseudo-elliptical NFW) . Dark to light 
colors delimit the confidence levels (from 1-a to 4-cr). 
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Fig. 8. Solid lines: distribution of the best-fit velocity dis- 
persion Vc from the optimisation of the lens configuration 
shown in Fig. ^ for each cosmological model. The cross 
(-I-) represents the original value for (ri°„, ri°) = (0.3, 0.7): 
Vc = 2000 km s"'^. Dashed lines correspond to the Icr con- 
fidence level contours from Fig. 0. 



3.3.2. Changing the shape of the mass profile 

To test the sensitivity of the method to the chosen fiducial 
mass profile, we tried to recover the lens with another po- 
tential, namely an elliptical HK profile, keeping the same 
simulated lens, tio, Oa and 6s were left free for the optimiza- 
tion. We first optimized the geometrical parameters for an 
arbitrary choice of cosmological parameters. The best val- 
ues found are: = 0.059", yo = 0.063", 6 = -0.063°, 
and e = 0.280. These values are close to the generating 
ones {xq — uo = 0", 6 = 0°), except for the ellipticity 
which does not correspond to the same physical mean- 



ing in the pseudo-elliptical NFW profile ( Golse fc Kneit 



2002). They were then kept fixed for the rest of the opti- 



1230 



-1-50 
50 



I90"tlo- The confidence 



mization. For the lens parameters, we found ctq 
km s^^ 

levels in the {ilm, ^\) plane are displayed in Fig. ^ 

Although the reconstruction with a potential model 
different from the initial "real" one does not perfectly 
fit the data, the results are quite satisfactory. The con- 
fidence levels are even tighter than in the previous case, 
but the HK-type potential is characterised by one addi- 
tional parameter or equivalently one degree of freedom less 
(i/ = 11), compared to the pseudo-elliptical NFW profile. 
Nevertheless we find a minimum reduced x^ — 5 rather 
far from 0. 

Several other mass profiles were tested as we wanted 
to discriminate between the different families of density 
profiles and test their sensitivity in the estimate of the 
cosmological parameters after the lens reconstruction. We 
used 5 types of profiles, namely: 

i) the pseudo-elhptical NFW profile ( poise fc Kncib 2002| ), 



ii) the singular isothermal ellipsoid (SIE) with p{R) = 
Pq/B?, R being the elliptical coordinate (Eq.(pl[)), 

iii) the isothermal ellipsoid with core radius (GIF), ob- 
tained by replacing R by R^ + in the previous ex- 
pression (see Kovncr 198_9| ), 

iv) the HK profile (Eq.(|)), 

v) and the King profile characterised by 



p{R) = pq 



R^/a" 



(l + i?2/a2) 



2+a ■ 



(26) 



The first 2 profiles are cusped, while the latter have 
a core radius and then an additional parameter. For each 
mass model, we generated the system of images defined in 
Table ^ (except for the SIE for which the radial system 
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Table 3. Results of the lens reconstruction using a mass 
model different from the one used to generate the systems 
of images. The minimum reduced is given for each sim- 
ulation. 



Input profile HK King CIE NFW SIE 

Fitted profile 
HK {u = 11) 
King (i/ = 11) 
CIE (f = 11) 
NFW {u = 12) 
SIE (i/ = 12) 



Table 4. Recovering of the free parameters of the lens 
potential for the Table || different systems of images. The 
errors represent the variation of the fitted parameters at 
1-tT level when scanning the {^m,^\) plane in the opti- 
misation process. 
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consists only of 2 images). We then fitted these images 
with the other 4 models. All the lens parameters were 
left free in this optimisation to get the minimum reduced 
X^- We did not change the cosmological parameters in 
these recoveries. The results are presented in Table ||. We 
note that the "core-radius" profiles (especially the HK and 
King ones) can easlily recover the systems generated by 
any other models. Indeed in the fit of cusped lens images 
by shallower profiles, the core radius can be reduced to 
very small values to mimic a large density slope near the 
center. This is not the case for the cusped models which 
cannot mimic images given by a finite core radius lens 
model. 



3.3.3. Influence of the number of multiple systems 

In the preceding sections we considered 3 systems of multi- 
ple images. As the method proposed is based on the differ- 
ence of angular distance ratios for different rcdshift planes, 
we now investigate the influence of the number of image 
families. The potential model is again an HK-type profile 
at zl = 0.3 with (To = 1400 km s-\ Oa = 13.5"(i.e. 65 
kpc), 0s = 146 "(i.e. 700 kpc) and e = 0.2. With 2 sys- 
tems of images, we consider only 2 free parameters for the 
cluster, because there are not enough observables to yield 
results for more parameters, while in the other cases, 3 
parameters are fitted. In all cases, these parameters are 
strongly constrained by the fit. Table ^ reports the errors 
on the fitted parameters in the optimisation process, for 
the different sets of multiple images detailed in Table ||. 
The differences in the fitted parameters between the dif- 
ferent cases are small, as they are already well constrained 
with a single multiple images system. 

The expected constraints on tighten when 

the number of families of multiple images increases 
(Fig. p^ ), especially when their rcdshift distribution is 
wide. 2 families would only provide marginal information 
on the cosmological parameters whereas 4 spectroscopi- 
cally measured systems would give very tight error bars, 
provided they are well distributed in redshift. 



Table 5. Sets of multiple images used in the simulations to 
test the influence of their number, rii represents the num- 
ber of images used for each family. It does not include the 
central de-magnified image created for tangential images. 
Nc is the number of constraints in the lens modeling for 
each family. Nc = 2 x — 2. 
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rii 
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1. 
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Tangential 
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Radial 
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3.3.4. Influence of additional galaxy masses 

In the previous parts, we considered only a main cluster 
potential with a regular morphology. We now test the con- 
tribution of individual galaxies, following the prescription 
used by Natarajan fc Kncib (1996 ) as in Section |3.l| . We 
generated 3 systems of multiple images formed by the sum 
of a main cluster with the mass density (HK-type) charac- 
terised by (To = 1400 km s-\ Oa = 13" and 0^ = 150" and 
12 individual galaxies which represent 30% of the total 
cluster mass (Fig. pi] ). 

The images were reconstructed using a main cluster 
potential with the same kind of shape as the initial one 
and the contribution of the galaxies scaled with cr^. In 
addition, we fixed ctq^ proportional to cto to avoid an in- 



crease of the number of free parameters. Consequently, 
any variation in ctq means a rescaling of the total mass of 
the cluster. So at first order we find that ctq E is constant 
when we scan the {ftm, ^\) plane. Keeping the geometri- 
cal parameters fixed {xq — Vo = 0", = 0°, and e = 0.2), 
we obtain the confidence levels in the {i^m, ^\) plane plot- 
ted in Fig. ^ and the following constraints on the poten- 
tial parameters: ctq 



1400^65 km s- 



13"^o:3 and 



= 151 



ii+i 
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Fig. 10. confidence levels in the (ri,„, JIa) plane obtained from the optimisation of the lens configuration described 
in Table ||. Left: 2 systems and = 10 — 4 = 6 degrees of freedom. Middle: 3 systems and v — 16 — 5= 11 degrees 
of freedom. Right: 4 systems and v — 2^ — h — \h degrees of freedom. The cross (+) represents the original values 
= (0.3, 0.7). Dark to light colors delimit the confidence levels (from 1-cr to 4-cr). 




Fig. 11. Multiple images generated by a cluster at zl = 
0.3 with an elliptical HK lens profile and the parameters: 
(To = 1400 km/s, Ba = 13.54"(a = 65 kpc) and = 
145.8"(s = 700 kpc). 12 individual galaxies are added in 
the potential. 3 families of multiple images are identified 
(see Table || for details) . We represent the radial (inside) 
and tangential (outside) critical lines corresponding to the 
multiple images redshifts. Their characteristic radii are 
increasing with redshift. Units are given in arcseconds. 

To test the influence of the individual galaxies, we tried 
a reconstruction without their contribution. For the geo- 
metrical parameters first optimised we obtain xq = 0.227", 
j/o = 0.060", Q = -0.748° and e = 0.193, stiU close 
to the generating values. The confidence levels in the 
(ilm,J7A) plane are plotted in Fig. |l^. The contours are 
slightly shifted and widened compared to the "good" ones 
(Fig. ^) but not significantly different. The minimum 
reduced is 17. So we are able to correctly retrieve 
the cluster potential, even without the individual galaxies 
(o-Q 
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Fig. 12. confidence levels in the ($!,„, JIa) plane ob- 
tained from the optimisation of the lens configuration 
shown in Fig. 
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1380l^g km s- 



9a = 11.9"^o:2 and 9, = 180"+3). 



For the individual galaxies, we assumed 
that their mass is scaled with the total mass with ct^ oc cto- 
The 3 ma,in cluster pctranieters (Tq, Oa and 9s were recov- 
ered with Xmin — ^^'^ V = 11 degrees of freedom. Dark 
to light colors delimit the confidence levels (from 1-ct to 
4-ct). 



Adding their contribution is nevertheless useful to deter- 
mine precisely the minimum region and to tighten the 
confidence levels. It becomes quite critical in more com- 
plex cases or when a single galaxy strongly perturbs the 
location of an image. 

3.3.5. Bi-modal cluster mass distribution 

Up to this point, we have considered simple clusters, domi- 
nated by a single massive component. In reality, most clus- 
ters are not fully virialised and present sub-structure as 
the result of accretion processes or merging phases. With 
these more complex mass distributions, the lensing config- 
urations are more widely distributed. Therefore we exam- 
ine how the cosmological parameters can be constrained 
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Fig. 13. confidence levels in the {yirm^x) plane ob- Fig. 15. confidence levels in the {Vlm,^\) plane ob- 



tained from the optimisation of the lens configuration 
shown in Fig 11. Here, we did not introduce the individual 
galaxies when recovering the global potential. The 3 main 
cluster parameters ctq: and 9^ were recovered but with 
a non-zero reduced x^in (Xmin ^ I''')- Dark to light colors 
delimit the confidence levels (from 1-a to 4-ct). 



tained from the optimisation of the lens configuration 
shown in Fig |lj. The main cluster parameters am, ao2, 
dai and 6*02 were recovered with Xmin — ^^'^ ^ number 
of degrees of freedom i/ = 6. The cross represents the 
original values = (0.3,0.7). Dark to fight colors 

delimit the confidence levels (from 1-a to A-a). 




Fig. 14. Multiple images generated by a bimodal cluster 
at zl ~ 0.3 with the lens parameters: (Tqi = '^02 = 1100 
km/s, Oai = 9a2 = 12" (58 kpc) and 0si = 9^2 = 167" 
(800 kpc). 3 families of multiple images are identified at 
zsi = 0.7, zs2 = 1 and 2:53 = 2. Units are given in arcsec- 
onds. 



with this type of realistic mass distribution. We thus gen- 
erated a bi-modal cluster consisting of two clumps of equal 
mass and 3 families of multiple images probing each part 
of the lens (Fig. |l^). The total potential is axisymmet- 
ric and each clump is characterised by an HK-type ellip- 
tical mass profile. As the number of multiple images is 
rather small, we limited the number of parameters to re- 
cover and chose ctq and 6a for each clump as adjustable 
variables. Therefore we fixed xqi — —34", xq2 — 34", 
yoi = Vm = 0", Oi = -45°, 02 = +45°, ei - 63 - 0.2 
and 9 si = 9s2 = 167". 

Fixing again the initial values of (il^„,il^) to the 
ACDM model (0.3,0.7), we obtain the confidence levels 
plotted in Fig. Ha. The contours are widened compared to 



Fig. 16. Multiple images generated by a cluster at zl = 
0.3 consisting of a main clump ((Tq = 1400 km/s, 9a = 
13.54" - 65 kpc - and 9^ = 145.8" - 700 kpc) and a smaller 
one (cto = 500 km/s, 9a = 5.2" - 25 kpc - and 9s = 45.9" 
- 220 kpc) located 102" from the main one. Close to their 
respective critical lines, 3 families of multiple images are 
identified: a tangential one (# 1, zgi — 0.6), a radial one 
(# 2, Z52 = 1) and another tangential one (# 3, zsz = 2). 
Units are given in arcseconds. 



the case of a single potential (in this case, the number of 
degrees of freedom is reduced from 11 to 6, but they still 
give reasonable constraints). Moreover we note that there 
is little variation in the fitted parameters: <toi = HOOiso 



km s-i, = 1100+^^ km s 



= 12.1 



//+0.1 

-0.1' 



and 



12.1 



//+0.3 

-0.2- 



This configuration is close to the case 
of the cluster Abell 370, modeled with a bi-modal mass 



distribution (Kneib et al. 19931; Bezecourt ct al. 1999b) 
needed to reproduce the peculiar shape of the central 
multiple-image system. Unfortunately, up to now only two 
redshifts are known for the multiple images identified in 
A370! 
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Fig. 17. confidence levels in the {yirm^x) plane ob- 
tained from the optimisation of the lens configuration 
shown in Fig. The 3 main cluster parameters (To: and 
were recovered with a reduced Xmin = 9 and a number 
of degrees of freedom v = 11. The cross (+) which rep- 
resents the original values (f2^,r2°) — (0.3,0.7) is now 
outside the 3-(T confidence levels. Dark to light colors de- 
limit the confidence levels (from 1-a to A-a). 



Last, we generated another system of 3 families of mul- 
tiple images produced by a cluster consisting of a main 
clump ((To = 1400 km s"'^) and a smaller one ((Tq = 500 
km s^^) representing 22% of the total mass (Fig. |l^). We 
chose to miss the small clump in the mass recovery as 
this may happen when dealing with some "dark clumps" . 
Fitting the configuration with a single main cluster, we 
found in a first round the geometrical parameters, which 
then remain constant in the ^^-optimisation: xq = 0.348", 
j/o = 0.189", e = 1.880° and e = 0.259. We note in par- 
ticular that the ellipticity is larger than the one used to 
generate the main clump (e = 0.2). This seems to be the 
response of the fitting process in order to mimic the miss- 
ing second clump. 

The parameters left free are again ctq, 9a and Og- The 
confidence contours are shown in Fig. |l^. We found the 
following values of the parameters: ctq — WOOI^q km s~^, 
9a = 12.8"^o'2 and 9s = 169"^2- However in this case, we 
do not recover correctly the set of cosmological parameters 
used to generate the system: (^I^j^a) — (0.3,0.7) is ex- 
cluded at the 3-cr level. Moreover the shape of the contours 
is not the one expected from the lensing degeneracy. This 
could be considered to be a signature of an incorrect fidu- 
cial mass distribution due to a missing clump in the mass 
reconstruction. This example demonstrates that the initial 
guess and the modeling of the different components of a 
cluster are very sensitive elements. They need to be care- 
fully determined if one wants to test further constraints 
on the cosmological parameters 



Table 6. List of 6 redshift configurations used in the com- 
bination of different cluster-lenses (Fig. |l^) for a global 
minimisation. 
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ZS3 


0.15 


0.4 


0.8 


2.0 


0.2 


0.5 


1.0 


3.0 


0.25 


0.6 


0.9 


2.0 


0.3 


0.6 


1.0 


2.0 


0.35 


0.6 


1.5 


3.0 


0.4 


0.8 


1.8 


4.0 



4. Conclusion and future prospects 

In this paper we have explored in detail a method to obtain 
informations on the geometry of the Universe with gravi- 
tational lensing. It follows an approach first presented by 
Link & Pierce (LP98) which states that multiple imaging 
systems at different redshifts can provide constraints not 
only on the mass profile of the lensing cluster but also 
on second order parameters like f2„i or Q\ - contained 
in angular size distances ratios. We have shown that this 
technique gives constraints which are degenerate in the 
(ri„j, plane and that the degeneracy is roughly perpen- 
dicular to the degeneracy issued from high-redshift super- 
novae searches. Moreover, the matter density can be 
better constrained than the A-term. Several simulations 
of lensing configurations are proposed, assuming reason- 
able conditions on the cluster-lens potential, such as a 
regular morphology modeled with only a few parameters. 
Provided high quality data can be obtained on at least 
3 systems of multiple images, such as high resolution im- 
ages (HST-type) for accurate image positions and deep 
spectroscopic data for the measurement of the source red- 
shifts, we can expect typical error bars of = 0.30±0.11, 
Vtx = 0.70±0.23. 

It is important to underline that one cluster-lens 
with adequate multiple images would provide by it- 
self a strong constraint on the geometry of the whole 
Universe. Such clusters are not that rare: MS2137. 3-2353, 
MS0440. 5-1-0204, A370, A1689, A2218, AC114 are cer- 
tainly good candidates for such an experiment. A thor- 
ough and detailed analysis is still to be done and we have 
in hand most of the tools to address the problem imme- 
diately. Furthermore, as the exact degeneracy in the (f2„j 
, S^a) plane depends only on the values of the different 
redshift planes involved, combining results from different 
cluster-lenses can tighten the error bars. For illustration, 
we combined 6 different lens configurations and source red- 
shifts, as listed in Table Compared to the expected re- 
sults with a single cluster (solid lines) , the constraints can 
be improved significantly (Fig. p^ . 

Looking for a good accuracy on the cosmological pa- 
rameters is a permanent search in cosmology. Although 
the curvature is now determined with a remarkable pre- 
cision thanks to recent results from CMB balloon exper- 
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Fig. 18. Color scale: x^{^rn,^x) confidence levels ob- 
tained for a combination of 6 different cluster-lenses con- 
figurations (see Table ^ for redshift informations). The 3 
main cluster parameters (Jq, 9s were recovered for each 
cluster with a reduced Xmin ~ ^'^'^ v = degrees of 
freedom. Dark to light colors delimit the confidence 
levels (from l-cr to 4-(t). Solid lines: confidence 
levels (from l-cr to 4-(t) obtained for a single cluster at 
zl — 0.3 (same as Fig. |^). The cross (-I-) represents the 
original values {yil^Ml) = (0.3,0.7). 



ime nts, it is still very diffic ult to disentangle from 
r^A ( Zaldarriaga et al. 1997 ). Therefore the advantages 
of joint analyses by se veral indepen dent approaches have 
been pointed out (see White (1998| ) and Efstathiou et al 
(1999|]]: combined results from the m — z relation for SNIa 



and CMB power spectrum analyses (which have orthog- 
onal degeneracies) constrain Q.m or Q.\ separately with 
much higher accuracy than the individual experiments 
alone, leading to the currently favored model. One impres- 
sive example has been given by Hu & Tegmark (1999) who 
showed that a relatively small weak lensing survey could 
dramatically improve the accuracy of the cosmological pa- 
rameters measured by future CMB missions. The combina- 
tion of independent tests can improve the constraints as 
well as serve as a consistency check. This is clearly demon- 
strated by Helbig et al. (1999) who combine constraints 
from lensing statistics and distant SNIa to get a narrow 
range of possible values for fi^. Therefore, gravitational 
lensing is a powerful complementary method to address 
the determination of the geometrical cosmological param- 
eters and probably one of the cheapest ones, compared 
to CMB experiments or SNIa searches. Our technique, 
when applied to about 10 clusters, should be included in 
such joint analysis, to obtain a consistent picture on the 
present cosmological parameters. We are truly entering 
an era of accurate cosmology, where the overlap between 
the allowed regions of parameter space is becoming quite 
reduced. 
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